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Since overhead costs for indirect labor account for a 
large percentage of the Naval Air Rework Facility's (NARF) 
total budget, it is essential that management be able to 
predict these costs accurately. The research performed in 
this thesis uses data from the major cost centers which 
comprise NARF Alameda. Regression models of their indirect 
labor to be used for forecasting purposes were developed. 
Quarterly data were used in the analysis, requiring 
transformation of the data to eliminate the effects of 
autocorrelation. The Durbin-Watson test was used to 
check for the effects of first-order autocorrelation and 
Wallis’ test was used for fourth-order autocorrelation. 

Once the effects of autocorrelation were eliminated, 
excellent structural results were obtained for twelve of the 
thirteen cost centers of interest. Predictive analysis 
performed using withheld data showed the final models can be 


expected to yield reliable forecasts. 
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ie JNERODUCTION 


The Naval Air Rework Facility (NARF) at NAS Alameda is a 
very large industrial rework facility which employs 
approximately 4500 civilians. These personnel are 
classified and budgeted as direct labor and indirect labor. 
The NARF operates on a fiscal budget which is evaluated and 
updated quarterly. It is imperative that personnel 
requirement forecasts be made accurately to ensure that the 
facility operates within its allotted budget. To this end, 
management is currently responsible for determining indirect 
labor requirements within their respective cost centers. A 
forecasting model to assist management in their decision 
making process would be of great value. The purpose of this 
thesis is to develop such a model. 

The objective of the author was to obtain data from NARF 
Alameda and analyze it with the goal of developing a number 
of forecasting models, each one unique to a major cost 
center, to be used by NARF management as a tool for 
determining future indirect labor requirements. The 
objective is a difficult one due to the complexity of the 
NARF organization and the necessity for parsimony in the 
models as they will be used by personnel from varying 


backgrounds. The models must be understandable, believable 


by management, and certainly they must be shown to be 
reliable. 

The development of these models requires a thorough 
understanding of the NARF organization, past and present, to 
ensure data are properly adjusted to reflect the 
Organization as it exists today. Understanding the history 
of the NARF command's interest in these models is also 
useful. An existing prediction model that is no longer 
being used iS examined to help in developing an 
understanding of the underlying feelings of the NARF 
management, as they are the group who will ultimately decide 
to what extent the models are utilized. 

Autocorrelation is almost always found in seasonal 
data. The analysis performed in the development of the 
indirect labor prediction models deals directly with the 
autoregressive process. Each model is used to predict 
indirect labor requirements for the last four quarters of 
available data. These predicted values are compared to the 
actual indirect hours worked. The method used for this 
predictive analysis is discussed. Results are presented 


with an explanation of the reliability of each final model. 
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If. BACKGROUND 


A. NARF ORGANIZATION 

The Naval Air Rework Facility (NARF) at NAS Alameda is 
an industrial activity of the Naval Shore Establishment that 
is under the command of the Naval Aviation Logistics Center 
(NALC). The NARF is a rework facility whose work includes 
repair, overhaul, conversion, modernization, Standard Depot 
Level Maintenance (SDLM), and analytical rework on 
designated weapons such as aircraft, missiles, engines, 
components, and associated accessories and equipment. NARF 
Alameda is a large organization which employs about 4500 
Civilians. Funds for operations (including overhead costs), 
manpower ceilings, equipment and tooling, and material 
Support are controlled and provided to the NARF by the 
Commander, NALC. Among other things the NARF commanding 
officer is accountable for the economy of operations, local 
management adaptations and adjustments, and maintaining and 
improving management. Periodic rework of aircraft and other 
weapons is scheduled into the NARF by NALC based on the 
Naval Air System Command's and Chief of Naval Operations’ 
calculated long range requirements. (NARF ALAMEDA INST. 


peo ia4d, 1974) 


ell 


The first element of Organization at the NARF is the 
command level, followed by the officer billets at the 
directorate level above departments. The next level 
includes departments which are subdivided into divisions, 
branches, sections, and units for service or production 
departments. A five-digit code identifies each of these 
work segments. The codes are commonly reduced to only the 
digits necessary for correct identification. Table 1 shows 
the codes and titles for the Command Office, Special 
Assistants, and Top Management. Table 2 shows an example of 
how the other management deSignations utilize the five digit 
codes to enable the exact identification of the 


organizational entity involved. 


B. COST ACCOUNTING PROCEDURES 

Volume I of the Cost Control Manual (NARF INST 7650.1) 
contains all the pertinent policies and procedures which 
govern the operation of the NARF's job order and cost 
control systems. It outlines direct and indirect 
expenditure guidelines, work classifications, and time 
keeping procedures for all NARF personnel. The codes used 
for the job order system are closely related to the 
Organization codes. The job order system is based ona 
concept of associating costs with end-products and overhead 
functions. Job orders are categorized as either direct or 


indirect. 
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OmEeucial Coding 


TABLE 1 


TITLES AND CODES FOR NARF ORGANIZATION 


Command Office 
00000 Commanding Officer 


01000 


Executive Officer 


Special Assistants 


00100 
00200 
OAais0'0 
01200 
01300 
01400 
01500 
01600 
Ua 0.0 


Occupational Safety and Health Manager 
Aviation Safety Officer 

Office of Counsel 

Civilian Personnel Director 

Position Management Office 

Security Office 

Public Atfairs Office 

Resources Management 

Demo Project 


Top Management (Directorates and Departments) 


02000 
03000 
04000 
05000 
07000 
20000 
30000 


40000 
50000 
60000 
70000 
80000 
90000 


Production Officer 

Resources Management Director 

Quality and Reliability Assurance Officer 
Weapon System Manager (WSM) 

NAVAIR Engineering Support Officer (NESO) 
Management Controls Department Head 

NAVAIR Engineering Support Office (NESO) 

Dept. Head 

Quality and Reliability Assurance Department Head 
Production Planning and Control Department Head 
Production Engineering Department Head 

Material Management Department Head 

Flight Check Department Head 

Production Department Head 


Abbreviated Coding 


00 
OL 


O01 
002 
aie 
O12 
013 
014 
015 
016 
017 
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TABLE 2 


SERVICE DEPARTMENT (56221) AND PRODUCTION 
DEPARTMENT (94111) CODES 


5 60 2 2a 
Department - Production Planning and Control Dept. [ 
Division - Power Plant Planning Division---------- 
Branca - Engine Planning Branch------------------- 
Section - Status and Control Section----------------- 
Uneage = Process Control Cen €@ 0 > 93> pg ee ee 

Sede 
Department - Production Dept.--------------------- 5 
Division - Avionics Division---------------------- 
Branch - Accessories Branch----------------------- 
Section - Instrument Section------------------------- 
Unit - Instrument Shop------------------------------ 


Direct job orders are those which identify costs to an 
end product. Indirect job orders are those which accumulate 
costs that cannot be identified with, or are not readily 
aSSignable to, an end product. Similarly there is direct 
labor cost and indirect labor cost. The Cost Manual (NARF 
INST 7650.1) defines direct labor cost as the actual payroll 
cost of the time spent by an individual which can be 
identified with an end product or service. An indirect 
labor cost iS a general expense (overhead) incurred by 
various cost centers, service departments, and other costs 
associated with services received on a plant wide basis such 
as supply, data processing public works, etc. The 
overwhelming portion of the direct labor costs come from the 


production work centers in the 900 department and most 
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indirect labor comes from the remaining cost centers 
(Service departments). The cost centers for NARF Alameda 
are shown in Table 3 and are referred to by their 


abbreviated codes throughout the thesis. 


TABLE 3 


NARF ALAMEDA COST CENTERS 


Organization Code Abbreviated Code 
00000 Command and Staff 007.01 
01200 Civilian Personnel Office 012/CPD 
05000 Weapon Systems Manager 050 
20000 Management Controls Department 200 
30000 NAVAIR Engineering Support Office 300 
40000 Quality and Reliability Assurance Department 400 
50000 Production Planning and Control Department 500 
60000 Production Engineering Department 600 
65000 Plant Services Division 650 
70000 Material Management Department 700 
80000 Flight Check Department 800 
90000 Production Department 900 
90300 Production Training Branch 903 
25000 NARF General 250 
93000 Metal and Process Division 930 
94000 Avionics Division 940 
95000 Airframes Division 950 
96000 Power Plant Division 960 


—_— oe ame eee cme cee cme cee cee qe cee ce es es es ee es es es es es ee es se ee ee es es es es es es es es es es es ee es ii ee 


Direct labor at the NARF is considered to be any 
position where the employee spends 50 per cent or more of 
his work hours over a year's period on a job which has a 
direct labor cost (direct job number). An employee's 
poSition is indirect if it is not direct. Indirect labor is 


further divided into two subcategories: fixed indirect and 
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variable indirect. Variable indirect labor is defined as an 
employee that is not direct but whose position supports 
direct labor and whose position varies somewhat as the 
direct workload varies. Fixed indirect positions are 
relatively unaffected by changes in direct workload and not 
immediately supportive of direct work. Management 
ultimately determines whether a position is direct or 
indirect, fixed or variable. 

The NARF's Funding Budget is the annual workload plan 
which identifies man-hour and dollar requirements and 
quantity of aircraft, engines, and other products to be 
inducted during each quarter of the budget year. It also 
provides the necessary information for developing the 
Operating Budget which is used for planning the day-to-day 
operations and sets forth goals against which actual 


performance may be measured. 


C. STATEMENT OF THE PROBEEM 

NARF Alameda attempts to work within the financial 
constraints of the Funding Budget. This workload prediction 
is approved by higher authority (NALC) based on the assigned 
direct workload. The amount of indirect labor needed to 
Support this direct workload is predicted by the NARF and 
shown as a ratio of direct to indirect labor levels. It is 
incumbent upon management to ensure that this ratio, which 


is based on predicted indirect needs, be met or improved 
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upon by reducing indirect labor over time and by accurately 
DEedietingaquarterly indirect hours in SrdermeEems tay within 
the Funding Budget. 

To date there has not been a reliable, proven way to 
predict required indirect labor for any time period in the 
future. The need to accurately predict indirect labor 
levels is recognized by the command but the solution is not 
Saovyemeect is difficult €O even approach the problem due to 
Gillessize of the organization, its accounting procedures, and 
its method of keeping records. Indirect costs account for 
approximately half of the NARF's total budget. 

The NARF employs about 4500 civilians at any given 
time. In an organization of that size there are many 
divisions and shops or units within each department. As 
time goes on there is development of new technology, change 
in management, and changes in assigned workloads, each of 
which contributes to the necessity of restructuring the 
Organization. These changes are often small relative to the 
total workforce and occur with little or no documentation. 
At NARF Alameda, many changes were made to the lower levels 
of the organization that are not possible to trace by 
looking at old memos, documents, charts, phone directories, 
or by talking to the "old timers". Even some of the major 
changes such as one which involved large external hiring and 
dissolution of an entire department are not possible to 


reconstruct in great detail. For example, in the third 
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quarter of 1983 the 100 department was dissolved and the 
012/CPD office was simultaneously created. This affected a 
lot more than just those two cost centers but it was 
impossible to account for where each direct and indirect 
person involved went or came from or to even tell exactly 
how many personnel were affected by the change. This 
results in a tremendous problem when trying to use 
historical data for making predictions. 

The problem of predicting future required indirect labor 
was also made very difficult by the methods used to document 
indirect expenses, direct expenses, and leave costs 
throughout the NARF. Records of direct and indirect 
personnel by head count are incomplete and not detailed. 

The data is kept rather casually on a handwritten form. The 
record contains only the total numbers of personnel on board 
at the department level and does not distinguish between 
direct and indirect personnel. Although the information 
evaluated for this thesis contained data back to 1979, some 
of the quarters show only a NARF total for personnel on 
board that quarter. 

Records of direct, indirect, and leave hours paid are 
available on microfiche as far back as first quarter 1979 
but the record is a cost accounting statement and therefore 
contains negative hours. These negative hours are very hard 


to trace, making this source extremely difficult to use, 
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even for a trained analyst; but this was determined to be 
the only reasonable source of information on labor hours. 
The need for a model to predict indirect labor was 
discussed with top management at NARF. All agreed that a 
model was needed but confidence in the possibility of 
developing a useful model was mixed. Those who expressed 
doubt felt that a previous model, called the Personnel 
Budget Model (PBM), was aS good as a prediction model could 
be, but that it could not be trusted. They wanted a model 
that they could have more control over, that was more 
understandable, easier to use, and reliable. Members of 
management that did have confidence in the old model did not 
know why it was abandoned and suggested updating it and 
reinstating its use. The data base on which the old model 


was built is no longer available. 


Diz PREVIOUS AND CURRENT MODELS 

The total forecasted direct labor hours for the fiscal 
year are currently split among the cost centers down to the 
shop/unit level according to an existing model implemented 
on the NARF's mainframe computer. The program is called the 
Computerized Workload Planning and Budgeting System 
(CWPABS). There does not currently exist at NARF Alameda a 
reliable model for predicting future indirect labor 


requirements. 


eg 


In 1974 the NARF identified the need for a reliable 
method of predicting the amount of indirect labor required 
in any given quarter in the future. This problem was 
addressed in late 1975 by Commander Oleson, USN, Management 
Services Officer and Comptroller, in a memorandum to the 
NARF Commanding Officer (Oleson, 1975). 

Oleson (1975) developed a model for predicting the 
indirect requirements based on the direct work loaded in 
each cost center by the CWPABS. Oleson called the model the 
Personnel Budget Model (PBM). The PBM was used in its 
Original form for a short time. The model was viewed by the 
top management as difficult to understand and they began to 
not trust it. Some adjustments were made to the model by 
NARF employees before it was totally abandoned, although the 
modified PBM still exists today on their mainframe. 

It was intended to use the PBM for determining the 
personnel budget for each cost center, outputting enough 
detail to enable the managers to understand the effect of 
factors such as increased training. The model generated 
personnel requirements based on the direct workload in each 
cost center as determined by the CWPABS, and indirect ratios 
associated with the direct. It used a combination of direct 
work in the production divisions (930, 940, 950, 960), 
direct work within the cost center being evaluated, overtime 
ratios, leave ratios, and time required for cleanup, 


training, etc. The addition of the fixed, functionally 
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related people in the cost center fee eee eee er eentect 
indirect budget required for the cost center. Leave ratios 
were not determined by the model or by Oleson but rather 
were ratios that were (and still are) determined within the 
NARF by the 500 department uSing historical data. Oleson 
Said the remaining required ratios could not possibly be 
derived by the use of regression techniques over time, that 
a "vintage" period had to be determined, and that the ratios 
be calculated from the data of that period. This idea and 
the methods used to determine the "vintage" period resulted 
in a period of only one year's time, 1974, being used as the 
data base to calculate all the remaining ratios used in the 
PBM model. (Oleson, 1975) 

It was determined that the author would not try to 
modify or update the PBM in any way but rather would 
determine a new data base and develop an entirely new model 
which would meet the thesis objectives of determining a 
model that would be understandable by the intended audience 
(NARF management and user personnel) and would be trusted 
and easy to use in its final form. 

Since the PBM was abandoned, the problem of predicting 
the required indirect work force just one quarter in advance 
has become a very large burden on the NARF top management 
group. A good model is necessary to help management 
effectively meet their budget goals and requirements. 


Currently, indirect requirements are being predicted each 
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quarter for the successive quarter by the NARF Management 
Group (NMG), a group of the senior managers who are the 
department heads for the major cost centers which comprise 
the total NARF budget. Predictions are made ona 
combination of each manager's "needs", the total number of 
indirect personnel presently on board, and the prediction of 
required indirect personnel based on a model that exists as 
an oversimplified version of the old PBM. This "prediction" 
method is actually just a way of allocating total indirect 
personnel presently on board. About his PBM model Oleson 
stated, "a ‘shakedown' period is recommended to ensure that 
it represents reality to the fullest extent possible. To 
this end efforts will be made to obtain full concurrence 
from all cost center managers" (Oleson, 1975). The "model" 
that is now being used does not comply with this. 

The development of this current model began in 1984 by 
the NARF department heads (Memorandum for the Executive 
Officer, March 1984). The purpose of the meeting that 
resulted in the model was to establish a method for 
allocating indirect staffing, to develop a plan for 
achieving a 54/46 direct/indirect ratio, and to establish 
indirect staffing levels by cost center for the subsequent 
quarter. The model developed does indeed allocate the 
indirect staffing but it does not predict required 
staffing. Each department head established the indirect 


levels based on the current quarter fiscal year ceiling plus 
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each one's anticipated needs. These needs were discussed 
and once all the department heads woraed on their needs a 
percentage of cost center to total variable indirect labor 
was developed for each cost center. No justification was 
given for these percentages or for the agreed upon variable 
levels, yet it was agreed that these percentages would be 
used for establishing future variable indirect levels. 
These percentages are still being used. Fixed indirect 
labor is determined and identified by each department head 
for the respective cost center. This is not a prediction 
model but rather a model for allocation of on board labor 
for the subsequent quarter. There is no model presently in 
use at NARF Alameda for the prediction or allocation of 
indirect labor for more than one quarter in the future. 
Managers and supervisors (civilian and military) tend to 
use the size of the budget for which they are responsible as 
a measure of their performance. It is therefore only 
matural that the overall indirect budget grows over time 
Since management has control and responsibility for the size 
of the indirect labor force in relation to the direct 
workload. As an individual is promoted to supervisor, he 
feels he needs to hire a secretary. AS a Supervisor is 
promoted to a higher level he needs two supervisors and 
another secretary to replace him, etc. Without a good model 
it is impossible to accurately predict a future indirect 


budget. 
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It iS imperative that the model be easy to use and 
understandable. Since managers are held accountable for 
their decisions they must be able to understand the 
forecasting techniques underlying the model if they are to 
base their decisions on the model's predictions. It is 
crucial that the management understand the forecasting 
system (the model) being used. (Bowerman and O'Connell, 


19797 Deze 
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ee eOuLREMENTS 

A data base was required which was at least quarterly 
1tself, to be used for the development of an indirect labor 
prediction model. It was determined that data for indirect 
labor, direct labor, and personnel leave from certain cost 
centers was needed in order to build this model. The data 
set needed to be large enough to obtain predictions which 
could be tested for reliability. This special type of 
quarterly data comes from a process called a time series. 
Time series regression was chosen because it has the 
potential to not only explain the past, but to predict the 
future behavior of the variables of interest (Ostrom, 1978, 
Brose 

Data was required for all the cost centers shown in 
Table 2. Those are all the major cost centers and are the 
ones that management is currently interested in making 


indirect labor requirement predictions for. 


B. SOURCES 
The dependent variables of interest are indirect hours 
and leave hours within the major cost centers shown in Table 


2. Direct hours is used as the independent variable. The 
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best source of data appears to be historical accounting 
data. This data was on microfiche and covered all cost 
centers' direct, indirect, and leave hours for each quarter 
from the first quarter of 1979 to the last quarter of 1986, 
inclusive. This accounting data is called the 7310-68 Cost 
Center Statement. The NARF no longer receives the 7310-68 
On microfiche but the equivalent data iS available ona 
computer printed hard copy which it does receive. 

The 7310-68 statement contains data for the hours that 
were charged to each cost center at the NARF. The cost 
center statement is subtotalled down to the shop and 
division levels and totalled at the major cost center 
level. This made it very easy to extract the data at the 


desired levels for the cost centers of interest. 


C. ADJUSTMENTS 

The data in the Cost Center Statement was mostly 
complete with quarters being fiscal beginning October first 
each calendar year. The portions of the statement that were 
misSing or unreadable were extrapolated from data that was 
readable. The microfiche statement was new in 1979 and was 
not quite as detailed in early years as in subsequent years 
due to a breaking-in period for personnel making required 
entries. This resulted in dropping four data points for 
FY79 in the cases where the entries were zero or missing 


(050 and 906) 
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There were some disadvantages of using the 7310-68 
statement that are inherent in an SOSSMC ING statement. The 
major problem is that such statements contain negative 
hours. The hours are adjustments made when it is discovered 
that shop hours were incorrectly submitted and therefore 
incorrectly charged. These adjustments were not necessarily 
made in successive quarters and often not even within the 
Same cost center. For example the 94212 shop may have 
bicemrectly charged 1200 hours of indirect labor in™the 
first quarter of 1981 and realized in the third quarter of 
1981 that the hours should have actually been charged as 
1200 direct hours in shop 93234, a different cost center and 
direct hours instead of indirect hours. The actual 
adjustment might not have been made until the last quarter 
of 1981 with no indication on the microfiche as to what the 
negative hours mean. A lot of time was spent adjusting the 
data for the negative charges in all the cost centers. It 
was necessary to make all identifiable adjustments to 
accurately reflect actual hours worked, rather than hours 
charged, within each cost center in each quarter. 

Another problem was that there is a cost center (250) 
that does not physically exist and is there mainly to 
capture charges for unexplained hours worked, and to capture 
unidentifiable negative adjustments. It was not possible to 
adjust these hours in any way. This cost center also shows 


hours for the Youth Employment Program (YEP) which is funded 
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differently than the rest of the hours in the Funding 
Budget. YEP does not affect the RECuece prediction model. 
The hours in 250 in the microfiche are included in the total 
hours for 200 cost center; these were subtracted out so they 
would not influence actual hours worked in 200. Prior to 
the creation of 250 in 1983, cost center 700 was a dummy 
location that was used as 250 is today. When 250 began, 700 
became a physical cost center. It was not possible to 
adjust the data to reflect where the people that work in 700 
came from. 

Leave hours were also a problem in the development of 
the model. The cost accounting statement does not 
separately list direct leave hours and indirect leave 
hours. Indirect leave hours are assumed by the 500 
department at the NARF to be proportional to the indirect 
hours worked and to the amount of direct hours worked. 
Quarterly leave hours are predicted for each cost center by 
the 500 department. The accuracy of their predictions was 
not investigated in this thesis, although they are to be 
added to the indirect hours forecast by the models presented 
herein to obtain total indirect labor required. 

Once the data were corrected to reflect actual hours 
worked within each cost center an adjustment had to be made 
for organizational changes and periods of large hires. This 
was done to the maximum extent possible within the 


guidelines of only attempting to adjust for the major 
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identifiable organizational changes. There were two major 
adjustments that were made in this attempt to standardize 
the data to reflect the NARF organization as it exists 
today. The first adjustment was made to reflect the 
deletion of the 42610, 42330, and 66300 shops which moved to 
96248, 93136, and 94300 respectively in FY84. The second 
Major adjustment was a result of the dissolution of the 
10000 cost center and the creation of the 01200 cost center 
in FY85. Making adjustments for the changes in indirect 
labor resulting from this was very difficult but was 
attempted so 012 would have a data base longer than only 
eight quarters. To make matters worse not all the people 
went from 100 to 012. Some went to 200 and some went to 
600. Memoranda were researched and the data was adjusted as 
accurately as possible. There was also a lot of outside 
hiring at that time which was accounted for in the final 
adjusted data set. 

In FY85 the cost codes were redefined although many 
Stayed the same. This did not affect the direct vs. 
indirect hours, so no adjustments were necessary for this 
change. Cost center 903 was created in the third quarter of 
FY83 to capture indirect costs for production apprentice 


hours. The 903 hours were subtracted from 900 totals. 
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D. PLOTS OF THE DATA 

Time series plots of the sauweted data were performed on 
all the cost centers management is currently interested in 
making indirect labor predictions for. These plots were 
examined for outliers. The automatic rejection of outliers 
1s not always prudent. Outliers should only be rejected if 
they can be traced to particular causes such as errors made 
when recording the data (Draper and Smith, 1981, p.153). 
All apparent outliers were researched and discussed with 
NARF personnel in an attempt to determine their validity. 
Figure 3.1 contains the plots of the data from cost center 
400 which are typical plots. The final adjusted plots for 
all the data are shown in the Appendix A. The plots show 
adjusted straight time hours for direct, indirect and leave 
hours. The leave hours are total direct and indirect. As 
previously discussed, it was not possible to separate these 
types of leave hours. The leave hour plots show the strong 
seasonality of the leave data. Only Straight time hours are 
of interest because overtime hours are worked by the same 
people who work the straight time hours. The prediction of 
overtime hours is not of interest in forecasting the number 
of indirect people required to support the direct workload. 

All the plots of the data show a strong seasonal effect 
within years. There is a yearly cycle in which the 
quarterly hours of successive years fall in the same 


relative position with respect to the other three quarters 
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Figure 3.1 Time Series Plots of Adjusted Data 
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of their respective years. This 1S mOe surprrsingecome ne 
data is a time series of quarterly data. The seasonal 
fluctuation in the hours worked is due mostly to the effects 
of leave. It is easy to see that the first quarter 
observation in a cost center for direct or indirect hours is 
almost always lower than the quarter preceding it and 
following it. Since the NARF is traditionally closed during 
the Christmas season, the hours worked are lower in that 
quarter. The plots of leave hours supports this. 

The actual data for 012/CPD was limited to only eight 
quarters. The cost center was created in the third quarter 
of 1985 as a result of the 100 cost center being dissolved, 
external hiring, and the joining ne CPD. ‘The 012/CPD pilcGe 
of indirect hours in the Appendix A shows the data adjusted 
for these changes but it was not possible to make similar 
adjustments to direct or leave hours. Direct hours within 
00/01, 930, and 940 show a general decline in hours, whereas 
500, 600, 700, and 050 show a gradual increase over the 
available data. The 200 cost center's direct hours went 
through two cycles with a decrease in hours in the early 
80's, a long period of rising hours, and a new cycle 
beginning this year. There iS a Similar pattern in the 300 
direct hours. The trend in direct hours for 400 and 600 is 
the same with hours declining rapidly, staying fairly 


constant for a time, and now starting to increase again. 
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The variable of most interest is the dependent variable, 
indirect hours. There is a steady rise in indirect in 012, 
00/01, 600, and 700, but an overall decline in 200, 500, 
650, 930, and 940. The indirect hours in the other major 
BOstmaenters are Cyelicwe Cost) centers. 700, 02/800, and 903 
all have sixteen or fewer observations as they are more 
recently created. 

As previously noted, leave hours are extremely 
seasonal. They show a general decline in most cost centers 
although the decrease in hours is very gradual and steady. 
Leave hours are increasing in 200, 600, and 050 but only 
sarge ly . 

The most noteworthy observation is the strong 
seasonality that all the data display which is expected of 
“historical quarterly data from such a large organization. 
This time series data is the subject of the following 
chapter which concerns development of the indirect labor 


prediction models for the NARF cost centers. 
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IV. MODELING OF INDIRECT LABOR COSTS 


A. AUTOCORRELATION 

An important application of time series models is the 
forecasting of future values from current and past values. 
The use of the observations at some present or past time, t, 
to forecast a value at some future time ttl can provide a 
basis for economic and business planning, production 
planning, inventory and production control, and the control 
and optimization of industrial processes. When a model is 
based on time series data there are inherent disturbances or 
noise in the process (Box and Jenkins, 1976, pp.1-3). That 
is, the disturbance occurring at one point of observation is 
correlated with another disturbance. When observations are 
made over time, the effect of the disturbance occurring at 
one period is likely to carry over into another period. 
Disturbances which occur in one time period are not likely 
to stop abruptly at the end of that time period but rather 
linger on for some time after the occurrence. While the 
effect of the disturbance lingers on, other disturbances 
occur, adding to the effects of lingering disturbances. The 
Shorter the time periods between the creation of these 
disturbances the greater the likelihood that the past noise 


will continue to add to the disturbance. So, the shorter 
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the periods of observations, the greater the chance of 
encountering these autoregressive disturbances. 

We examine the reSiduals from a regression model (the 
"stochastic" disturbance) to determine if these error terms 
are related. If the disturbance occurring at time t is 
related to the disturbance occurring at time (t-s), they are 
said to be autoregressive. The presence of autocorrelation 
is more likely to occur when dealing with quarterly data 
than with annual data which is a much longer time period 
(Kmenta, 1971, pp.269-270). The residuals of a statistical 
model applied to quarterly data can be expected to exhibit 
autocorrelation of some form. Most applications of time 
series regression analysis assumes that a first-order 
autoregressive process is generating disturbances (Ostrom, 
1978, p.24). This model is called an AR(1) process. Kmenta 
(1971, p.271) considers a model in which the error terms are 


generated by the AR(1) process of the form 


Y= Fes, & SC, (4.1) 
and 


+ VY. ’ C422} 


where €,is the error term from a regression model 
corresponding to the observation at time t and €4., is the 
error term at time t-1. The » term is a coefficient of 


correlation between the related error terms, €,, and €, of 
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lag 1. The last term, W% , is a normal and independently 
distributed random variable with a zero mean and constant 
variance ag, that is assumed to be independent of &,., . 

The assumption often made in regression that the error 
terms are independent identically distributed normal random 
variables with a mean of zero and constant variance o° is 
not true with autocorrelation. The terms are in fact 
related to previous error terms and are dependent upon the 
form of autocorrelation which exists. In fact, the variance 
of the error terms is cme) (ies ) when the autocorrelation 
present is AR(1) (Kmenta, 1971, pp.271-272). 

The properties of the least squares estimators of the 
regression coefficients are also affected when the error 
terms are autocorrelated. The most important of these 
properties are bias, variance, consistency, efficiency, and 
linearity. When the error terms are autocorrelated the 
least squares estimators of the regression coefficients are 
consistent but are no longer efficient (Klein, 1974, 
pp.55-87). The estimators are also unbiased but the 
Standard estimates of the variances of the coefficients are 
biased. This affects the significance levels of the t and F 
tests (Bennett, 1979, pp.245-248). A positive 
autocorrelation in the error terms results in the estimated 
variance being underestimated. This could cause a serious 
overestimation of the t-statistics and significance levels, 


leading to unwarranted confidence in the regression model. 
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It iS easy to see that autocorrelated errors can be a 
serious problem in evaluating the fit of the regression 
model (Johnston, 1984 , pp.310-313). It is imperative 
therefore to obtain not only unbiased estimates of the 
regression coefficients, but to have unbiased estimates of 
their standard errors as well. 

As indicated in the AR(1) model, the error terms in any 
given period are related to those one period prior. 
Successive disturbances are frequently positively correlated 
in time series data (Theil, 1978, p.302), resulting in the 
adverse effects mentioned above. In addition, when data are 
quarterly observations, a special form of the fourth-order 
autoregressive process may be present. This special 


fourth-order autoregressive process has the form 
eee teeny tae ; (4.3) 


where the % are independent normal random variables, and 
the error terms are correlated with the errors in the 
corresponding quarters of successive years (Wallis, 1972, 
pp.617-621). The variance of the error terms when this 
special fourth-order autocorrelation is present is 

o} /(1- 7, ) (Judge et. al., 1985, p.298). The simple 
fourth-order autoregressive model shown in equation (4.3), 
3s henceforth called the AR(4) model. The AR(4) process 


used herein assumes that the effects of the past three 


oy 


consecutive quarters are negligible compared EO thesebfect 
of the corresponding quarter of the proves year (Boger, 
19837 bone 

In general the time series plots of the independent 
variable total direct labor from the 930, 940, 950, and 960 
production cost centers shown in Appendix A, Suggest that 
this AR(4) model is appropriate for these centers. Time 
series plots of the dependent variable indirect labor within 
each cost center, also shown in Appendix A, suggest the 


Same. 


Bs THEORETICAL MODEL USED 
The theoretical model used for determining indirect 


labor in this analysis is of the form 


Nan X,A+é (4.4) 


and 


oe Pewee + Ve ’ Ct = ite Noa eee (4.5) 


where in general X, is a T x k matrix, A is a k x 1 vector, 
and iis either 1 for the AR(1) process, 4 for the AR(4) 
process, or a mixture of the two processes. It was 
determined that k would be 2 in all cases to obtain 
parSimonious models that satisfactorily explain the 
dependent variable. Xz is a column of ones, for a constant 


term, followed by a column of the data for the independent 
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variable, total direct labor as previously described. The 
Beare the dependent variable, indirect labor within a 
specific cost center. The VY of the error term shown in 
equation (4.5) are independent normally distributed random 
variables with mean zero and constant variance. 

Results of equations (4.4) and (4.5) cannot be 
immediately used because when autocorrelation exists the 
estimators of the regression coefficients are not efficient 
and their variances are biased. Both X and Y can be 
transformed to eliminate the effect of the autocorrelation 
Powe to SNOW a regresstongen Che transformed sx and ¥ 
values, called the Generalized Least Squares (GLS) 
procedure, yields results that are corrected for the problem 
of autocorrelation. Since » is not known, it must be 
estimated from the sample observations (Kmenta, 1971, 


pome2—=285). 


C. DETERMINATION OF AUTOCORRELATION ESTIMATOR 

To make the transformation of the original data to 
correct for the known presence of the AR(1) or the AR(4) 
process in the errors, estimates of » and if necessary /p, 
are needed. The transformation for an AR(1) process is 
discussed in Judge et.al. (1985, p.285) and in Johnston 
(1984, p.318). The estimate of » suggested by Johnston 


(1984) is of the form 


So 


yf) =] ==e0, Saar ; (456) 
where 
= 
> (e.-e,.,) 
a *=55--- ieee} oe Lae (427 } 
=e, 


and the e, are residuals from the fitted least squares 
regression shown in equation (4.4) (Durbin and Watson, 1950, 
pp.424=425) 

Several estimators are available for ». In this thesis 


we use an eStimator of the form 


“~ 


Py = il as QK 5 d, a (4.8) 
where 
a 
> (.-e,-4) 
t:S$ 
Gh Se (4.9) 
Sie < 
ta) 


and the e, are the same as in equation (4.7) (Wallis, 1972, 
pp.617-624). This estimator has been shown to perform well 


(Boger, 1987, p.5). 


D. TRANSFORMATION FOR AUTOCORRELATION 
Once the appropriate estimates are calculated, the data 
can be transformed using the estimates. The transformation 


used for the AR(1) and AR(4) processes are 
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ae ee) ae E~ le... i 10) 
and 
¢ 
= eee, c= ile, (4.17) 
&N : . 
where »”, i = 1 or 4, is the estimate from (4.6) or (4.8) as 


; *& 
required. The Zs are the transformed data from the 


wiemEece and direct labor values. 


Pee EROCEDURE 

Each cost center was evaluated using the following 
sequence of procedures. The general procedure was to first 
perform the ordinary least squares (OLS) regression with the 
dependent variable being the particular cost center's 
indirect labor. In all cases the independent variable was 
the summation of the direct labor from production cost 
centers 930, 940, 950, and 960. Other predictor variables 
were evaluated and parameters, in addition to the intercept 
term and above predictor, were tried. It was determined 
that there was no Significant improvement that could be 
found in any of the models by using any form other than the 
parsimonious one containing one predictor. After completing 
an OLS regression the residuals were analyzed and tested for 
the presence of AR(1), Wallis's special AR(4), or a 
combination of both. This was done by looking at a plot of 


the autocorrelation function of the residuals to get an 


4l 


overview of the type of autocorrelation present, if any. 
More formal testing was then performed on the residuals. 

The Durbin-Watson test was used to check for the 
presence of AR(1) in the regression model (Durbin and 
Watson, 1951, pp.159-175). The Wallis test has been shown 
to be a generalization of the Durbin-Watson test and was 
used to check for the presence of the AR(4) process in the 
regression model (Wallis, 1972, pp.621-624). For both the 
AR(1) and the AR(4) cases a two-sided test was performed, 
using the null hypothesis of zero autocorrelation against 
the alternative of non-zero autocorrelation with a 
Significance level “< = 0.10. 

Both of these tests share the awkward problem of an 
inconclusive region. When d, < d < d, , as culculated in 
(4.6) or (4.8), the test is inconclusive and it cannot be 
determined if autocorrelation is present. This problem is 
accentuated when the sample size is small (Johnston, 1984, 
p.316). In the present study, the largest sample size is 34 
and the smallest is only 13. This problem was dealt with by 
following the statistically conservative procedure of using 
the upper significance point, d, , as if it were the 
critical value. The actual lower significance point, d, , 
is ignored completely. This method was used for both the 
Durbin-Watson tests and the Wallis tests. This procedure is 
described in the literature as performing well in many 
Situations (Draper and Smith, 1981, p.167). The test 


procedure is now to use the following as the rejection 
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Sr Crecmnemetne twos ided sbest; wif deqd, or.d > 4 - d, , 
reject the null hypothesis at the 2 x levelPuesit1S easily 
seen that any point that previously would have fallen in the 
inconclusive region would now fall in the critical region 
and lead to the rejection of the null hypothesis. This 
procedure is also recommended by Johnston (1984). He states 
that it is more serious to accept the null hypothesis when 
autocorrelation is present than to incorrectly assume it to 
be absent. He also notes that when the regressors are 
slowly changing series, which many economic series are, the 
true critical value will be close to the upper bound 
(Johnston, 1984, p.316). This all gives credence to using 
the upper significance point in performing these tests. 
Since the Wallis test is only a slight modification of the 
Durbin-Watson test, this rule was applied to the Wallis test 
as well. 

The next step in the general procedure was to transform 
the data using the appropriate procedure depending on the 
form of autocorrelation found to be present. When both 
AR(1) and AR(4) were present the transformation was 
performed using whichever form was determined to be 
dominant. The regression coefficients were then reestimated 
using the transformed data of the dependent and independent 
variables to obtain a GLS solution. The residuals of this 
regression were tested for the presence of autocorrelation 


just as was done in the OLS regression. This procedure of 
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reestimating the model and checking for the presence of 
autocorrelation was repeated until model was obtained for 
which the residuals appeared to be free of any 
autoregressive process. 

Once a final model was obtained, the residuals from this 
model were checked to ensure the assumptions that they were 
independent, identically distributed, normal random 
variables having a mean of zero and a constant variance 
could not be rejected. 

The final step was to generate predictions to test the 
model's prediction capabilities. This was done in all cases 
with the exception of the 700 cost center which had too 
small of an original sample size. The last four 
observations were withheld from the data and the estimates 
were calculated as described above starting with an OLS 
regression and ending with a model using data that had been 
transformed as necessary. These predictions were then 
compared to the actual values of the indirect hours which 
were withheld. This method of withholding data has been a 
matter of discussion for some time. It has been argued that 
the model selection procedures described above amount to a 
considerable data mining. It thus seems wise to save some 
of the sample data for use in validating the resulting 
model. It is further argued that if the model passes its 
predictive test evaluation the four withheld observations 


should then be incorporated back into the data set to 
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reestimate the parameters with all available observations. 
If this model shows that its Perecdst 5 are sufficiently 
accurate, the four withheld observations are combined with 
the rest of the original data and the model parameters 
reestimated. This final model is then used to generate 
fOomecadstsmmmto the future (Ostrom, 1978, pp.58-59). Boger 
(1983, pp.33-40) suggests this procedure for predictive 
analysis and uses it by comparing the predicted values of 
the dependent variable with the observed values uSing a 
Pearson correlation coefficient, the root mean squared 
forecast error, and the mean absolute percentage error. The 
results obtained for the applications presented in the 
following chapter include these values, along with an 


explanation of their meaning. 
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A. GENERAL 

In this chapter the procedures described in Chapter Four 
were followed to obtain GLS regression models for each of 
thirteen cost centers. Separate models were obtained for 
each of the cost centers 012, 00/01, 200, 300, 400, ° 500, 
600, 650, 700, 930, 940, 9SQ7 and 960. Due co ene fare 
number of models and the repetitive procedures, only the 
case for the 400 cost center is presented in detail. The 
final results are presented for each of the other cost 
centers. In all cases the summation of direct labor hours 
from the 930, 940, 950, and 960 production cost centers was 
used as the independent variable and indirect labor hours 
within the particular cost center was the dependent 
variable. The computer programs used in the structural 
analysis are shown in Appendix B. All the adjusted data 
utilized in the analysis is presented in Appendix C so the 
models can be maintained as new quarterly data becomes 
available. 

Models were not obtained for cost centers 02/800, 900, 
903, or 050. It is easily seen from the time series plots 
of these cost centers in Appendix A that the creation of a 


prediction model is unwarranted due to the small number of 
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indirect personnel in each center and the relatively small 
Fluctuation of indirect hours worked in each quarter. 
Although not part of the structural analysis, the adjusted 


data for these cost centers is provided in Appendix C. 


ye PROCEDURE 

Table 4 presents the results of the procedures described 
above aS applied to the regression of indirect labor for the 
400 cost center (ID400) on the summation of total production 
direct labor (DIRECT). The best model was obtained using 
the last twenty-five available observations. This initial 
regression provided fair results if one ignores the 
autoregressive possibilities. The adjusted R-squared value 
is somewhat low but the F-statistic (not including the 
constant term) was well above its five percent critical 
value of 4.28. It must be remembered that both of these 
statistics were inflated due to the presence of 
autocorrelation. Even though the standard errors of the 
regression coefficients were biased downward due to the 
presence of the autocorrelation they were still large 
compared to the magnitude of the coefficients. 

The autocorrelation function of the residuals obtained 
from this initial OLS regression (Figure 5.1), with a large 
spike at lag one, strongly suggested the presence of first 
order autocorrelation. Upon formally testing the residuals 


for the presence of AR(1) using the Durbin-Watson test 
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TABLE 4 
RESULTS FOR COST CENTER 400 


Model: ID400 = a + b DIRECT 


Untransformed Data 


Standard Error of the Regression: 6638. 
Adjusted R-squared: oo 
F-statistic (degrees of freedom): 34.5 (1,23) 
Estimate of a: -12833. 
Standard Error: 14204. 
T-statistic -.90 
Estimate of b:;: 099 
Standard Error: 0017 
T-statistic 5.87 
Durbin-Watson Test Statistic: ~/0 
Wallis Test Statistic: 1.40 
Estimator for First Transformation (Af) 6493 
Estimator for Second Transformation (A, ) ©4565 


Transformed Data 


Standard Error of the Regression: 4379. 
Adjusted R-squared: OZ 
F-statistic: 108.4 
Estimate of a: -655. 
Standard Error: 1734. 
T-statistic -.3/ 
Estimate of b: ~087 
Standard Error: 008 
T-statistic 10.4 


Prediction Results 


Correlation Coefficient (CC) Between 054 
Actual and Predicted Values 

Root Mean Squared Error (RMSE) Divided by 05 
the Mean of the Actual Values 

Mean Absolute Percentage Error (MAPE) 4.06 


(in percent) 


Actual Values 70078 69305 62138 63621 
Predicted Values 64623 65649 62122 65639 
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of the Residuals for Cost Center 400. 


statistic (Durbin-Watson, 1951, pp.159-175), the null 
hypothesis that no first-order autocorrelation was present 
was Clearly rejected. The Wallis statistic for AR(4) was 
Pmicwgnrticant at this»point. 

Next the data were transformed using the calculated 
estimate of fg and GLS estimates of coefficients in the 
regression model were obtained. The residuals with this 
reestimated model were analyzed for the presence of AR(1), 


AR(4), or both types of autocorrelation. The Durbin-Watson 
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Statistic was no longer significant but now the Wallis test 
Statistic was significant indicating the presence of 
fourth-order autocorrelation. The data were transformed 
again, this time to eliminate the AR(4) process, and the 
regression coefficients were reestimated. Both the 
Durbin-Watson and the Wallis test statistics were clearly 
insignificant for residuals with this model, thus the null 
hypothesis of the presence of either the AR(1) process or of 
the AR(4) process in the residuals was rejected. The 
estimated autocorrelation of these residuals also showed 
reduction in the spike at lag one. Figure 5.2 shows this 
autocorrelation function of the residuals from the final 
model which indicates that no autoregressive process 
remained. 

The residuals were formally analyzed to determine if any 
of the assumptions required for the regression could be 
rejected. Tests were done to check the normality of the 
residuals and to test the constant variance (homogeneity of 
variance) of the residuals. 

The first test to check for the normality of the 
residuals was to generate an empirical cumulative 
distribution function (CDF). This compares the CDF of the 
residuals with that of the normal distribution. The 400 
cost center appeared to have a normal distribution with a 
mean of zero and a standard deviation of 4287. Next a 


probability plot of the residuals (a Q-Q plot) was 
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400 FINAL 


AUTOCORRELATION 





LAGS 


Figure 5.2 Autocorrelation Function of the 


Residuals for Cost Center 400 Final Model 


performed. This plotted the quantiles of the the residuals 
against the normal distribution. These plots are shown in 
Figure 5.3 and indicate a reasonable fit. The plots were 
also bounded by the ninety-five percent Kolmogorov-Smirnov 
(K-S) confidence boundaries shown by the dotted lines. The 
K-S test statistic was not significant (<= .90). 
Therefore, the null hypothesis of a normal distribution with 
a mean of zero and a standard deviation of 4287 could not be 


rejected. 
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CUMULATIVE DISTRIBUTION FUNCTION 
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Two additional tests were performed to check the 
assumption that the residuals were homoscedastistic. The 
residuals were plotted against the predicted dependent 
variable to check for any visual abnormalities or obvious 
patterns (Figure 5.4). The plot shows no discernible trends 


and the variance seems constant throughout. A statistical 


4000 
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FITTED VALUE 
Figure 5.4 Test for Homogeneity of Variance of 


the Residuals for Cost Center 400 Final Model 


test was next performed to further help validate this. The 
procedure as presented in Mood, et.al. (1974, p.438) tests 
the null hy,pcthesis that two sample sets have equal 
variance. The residuals were randomly divided into two sets 


and the null hypothesis H,: o* = os* was tested against the 
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= 
alternative hypothesis H : o, # oa %° , where g;* and a; 
were the variances of the two sets Of residuals. The test 


Statistic for this test is 


= 2 
(n,-D ECKL - XK, ) 


"GE OR = 


which has the F distribution with (n, - 1) and (n, - 1) 
degrees of freedom under H,. With 12 and 13 degrees of 
freedom the 400 cost center had a test statistic R of 
1.105. The null hypothesis was not rejected so the 
assumption of constant variance could not be rejected. 

All these tests could not reject the hypothesis that the 
residuals are normally distributed random variables with a 
mean of zero and a constant variance. It is noted that not 
rejecting these assumptions does not mean that they are 
necessarily correct but rather, that on the basis of the 
data observed, there is no reason to say the asSumptions are 
incorrect (Draper and Smith, 1981, p.142). 

The final model was a great improvement over the initial 
OLS regression solution. It indicated there was indeed a 
good deal more information about the indirect labor than 
initially seen in the OLS regression solution. The 
R-sgquared value was significantly higher indicating that the 


model contained a lot more information than just the mean of 
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the dependent variable. The t-statistic indicated that the 
intercept term was not significantly different from zero so 
it was not included in the model. The standard error of the 
remaining regression coefficient, the slope term, was 
relatively small compared to its coefficient. 

After transforming the data twice to eliminate the 
presence of autocorrelation from the residuals a GLS 
regression model was finally obtained which yielded 
excellent results. This final model was 

ID400 = .0878 DIRECT. 
It 1S important to recall that this model uses data that has 
been transformed twice using the appropriate autocorrelation 
estimators pA, and /y, which are listed in Table 4 along 
with the other results for the model. 

The final step was to test the predictive capabilities 
of the model as previously mentioned. The general procedure 
was to reestimate the model as described above but using 
only the first T-4 observations. Then four predicted values 
of the dependent variable were compared to the last four 
observed values which were withheld from the original data. 
These last four observations were predicted uSing the 


eguation 
¥ = AY, + (Xe > AX VA »  teT-Bye00,T, (542) 
where y, are the known or predicted values for the dependent 


a2 


variable and X4 are the conditional (always known) values 
for the independent variable. The value of i is either 1 or 
4 depending on whether the AR(1) or the AR(4) process was 
eliminated, and 2 and p. are the values obtained from the 
structural estimation based on the first (T-4) 

observations. (Boger, 1987, pp.10-11) If the model 
required transformations to eliminate the presence of both 
first-order and fourth-order autoregression the predictions 


were made uSing 


~~ 


oe A AA A A A > 
YUP, Yet Py Sonya fyi) | | ea eee rere, Ay Xeng YB (5.3) 


where the variables are described above. The computer 
programs used to make the predictions are shown in 
Appendix B. 

The reliability of the predictions was measured using 


the Pearson correlation coefficient (Theil, 1978, p.85) 
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CC = ee eee 22 2 2 2 ee ee 2 e222 2 2 eee eee ee eee eos ? (5 oan) 


ie = (4e- 7) Gz (F,-4,) 32 


the root mean square error divided by the mean of the actual 


M. 


last four observations (Ostrom, 1982, p.66) 


RMSE = -=4 73) ------------ ; (Soar 
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and the mean absolute percentage error 
4 se 
NU a a= SS SS Oo . (5.6) 


The Pearson correlation coefficient (CC) value of .541 
indicates that there is a tendency for the predicted values 
of indirect labor hours to follow the actual values. The 
ratio of tne root mean squared error to the mean of the four 
actual values to be predicted indicates a meaSure of the 
Size of the forecast errors. The root mean Squared error 
(RMSE) for the 490 cost center was 5.2 percent of this mean, 
which shows that the forecast errors were small relative to 
the actual values. The other measure of the size of the 
forecast errors is given by the mean absolute percentage 
error (MAPE). This measure for the 400 cost center 
indicates that the forecast errors were about 4.06 percent 
of the actual observed values (Boger, 1983, pp.34-36). In 
Summary, the results of the final GLS model for the 400 cost 
center can be expected to produce reliable, accurate, and 
acceptable predictions to be used as a tool by NARF 
management to assist them in their indirect personnel 
planning for four quarters into the future. 

This same general procedure was completed for all of the 
remaining models. The final results are tabulated along 
with a brief interpretation of their meaning and 


usefulness. There were only two differences in the 


7 


procedure for these remaining cases. The first was the 
selection of the sample size used in the regression 
analysis. Different sample sizes were used for each model 
and the one which yielded the best results was selected. 

The second difference in the procedure was the order in 
which the autoregressive processes were removed from the 
residuals. The AR processes were removed by order of 
Significance. The autocorrelation function of the residuals 
For each initial model was plotted and examined (Figure 
5.5a, 5.565, and 5.5c) to see which existing form of 
autocorrelation was dominate (closest to one). This was 
combined with the results of the Durbin-Watson and ‘Jallis 
tests to decide which AR process the data were to be 
adjusted for first. Subsequent order was decided by the 
results of the Durbin-Watson and Wallis tests performed on 
the residuals of the preceding model. The AR(1l) process was 
dominant in every case and the AR(4) process was adjusted 
for in the final regression of the models for the 012, 300, 


500, and 600 cost centers. 


C. THE REMAINING MODELS 
le Cost Center_Ql2 
Table 5 contains the results for the 012 cost 
center. Very poor results were obtained for the initial 
regression. Although the estimated autocorrelation function 


of the residuals of this model (Figure 5.5a) suggest that 
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Figure 5.5a Autocorrelation Functions of the Residuals. 
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Figure 5.5b Autocorrelation Functions of the Residuals. 
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Figure 5.5c Autocorrelation Functions of the Residuals. 
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TABLE 5 
RESULTS FOR COST CENTER 012 


Model: ID012 = a + b DIRECT 


Untransformed Data 


Standard Error of the Regression: ioe Ole 
Adjusted R-squared: -.034 
F-statistic (degrees of freedom): 036 (1,28) 
Estimate of a: 51126. 
Standard Error: 59109. 
T-statistic EO On 
Estimate of b: 001 
Standard Error: ~006 
T-statistic 19 
Durbin-Watson Test Statistic: 92 
Wallis Test Statistic: ~ ie 2 
Estimator for First Transformation (Py ) - 3880 
Estimator for Second Transformation (2, ) -so0U 


Transformed Data 


Standard Error of the Regression: ST. 
Adjusted R-squared: 2514 
F-statistic: 55,606 
Estimate of a: 12823. 
Standard Error: T9325. 
T-statistic 6.66 
Estimate of b: 028 
Standard Error: 2005 
T-statistic SG 


Prediction Results 


Correlation Coefficient (CC) Between -.08 
Actual and Predicted Values 

Root Mean Squared Error (RMSE) Divided by 55 
the Mean of the Actual Values 

Mean Absolute Percentage Error (MAPE) 12.9 


(in percent) 


Actual Values 53553 49303 48932 54049 
Predicted Values 46883 48520 41461 42023 


62 


the AR(1) process is the most significant form of 
autocorrelation, the best model reegosund to result from 
first adjusting the data to eliminate the presence of the 
AR(4) process and then the AR(1l) process, Examination of 
the residuals of the final model showed no indication of 
autocorrelation present and that none of the necesSary 
assumptions pertaining to the residuals appeared to be 
Violated. The final model is a great improvement over the 
results from the initial regression. Virtually all the 
indirect hours in the Civilian Personnel Department (012) 
are fixed indirect which by definition do not vary with the 
direct. The prediction results are actually quite good when 
this is considered, 

goose _ Center 9079) 

The results for the model of the Command and Staff 
cost center (00/0]) are shown in Table 6, The data was 
transformed to eliminate the presence of AR(1) and the model 
was then reestimated. Analysis of the residuals showed no 
apparent indication of autocorrelation present nor 
violations of the required assumptions. Like 012, there is 
no variable indirect that can usually be identified within 
the cost center so it is not surprising that the model is 
unsatisfactory. 

3. Cost_Center 200 
The final model for the Management Controls 


Department (200) is presented in Table 7. The results from 
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TABLE 6 
RESULTS FOR COST CENTER 00701] 


Model: IDQ0001 = a + b DIRECT 


Untransformed Data 


Standard Error of the Regression: 62213 
Adjusted R-squared: ©9013 
Fe-statistic (degrees of freedom): 1 a (1,28) 
Estimate of a: 30953. 
Standard Error: 11230. 
T-statistic 2.16 
Estimate of b: -.016 
Standard Error: 014 
T-statistic -1.18 
Durbin-Watson Test Statistic: ~ 406 
Wallis Test Statistic: ~ 423 
Estimator for the Transformation (2 ) .7970 


Transformed Data 


Standard Error of the Regression: 30S. 
Adjusted R-squared: eS 
Festatistic: ~137 
Estimate of a: 3957. 
Standard Error: 968. 
T-statistic 4.09 
Estimate of b: sC02 
Standard Error: - COD 
T-statistic » 369 


Prediction Results 


Correlation Coefficient (CC) Between -.97 
Actual and Predicted Values 

Root Mean Squared Error (RMSE) Divided by 41 
the Mean of the Actual Values 

Mean Absolute Percentage Error (MAPE) SJ cock 0 


(in percent) 


Actual Values 23433 25140 25585 27186 
Predicted Values 19554 16725 14264 12471 
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TABLE 7 
RESULTS OR. COST CENIER. 200 


Model: ID200 = a + b DIRECT 


Untransformed Data 


Standard Error of the Regression: 4103. 
Adjusted R-squared: 7 
F-statistic (degrees of freedom): eee ns (1,28) 
Estimate of a: 24720. 
Standard Error: 7406. 
T-statistic 3.34 
Estimate of b: 038 
Standard Error: 009 
T-statistic 4.25 
Durbin-Watson Test Statistic: 0994 
Wallis Test Statistic: 1.67 
Estimator for the Transformation (f/ ) a50e2 


Transformed Data 


Standard Error of the Regression: a5). 
Adjusted R-squared: 64 
F-statistic: 51.8 
Estimate of a: 10495. 
Standard Error: 2597. 
T-statistic 4.04 
Estimate of bs: 043 
Standard Error: ~006 
T-statistic oe 


Prediction Results 


Correlation Coefficient (CC) Between -.595 
Actual and Predicted Values 

Root Mean Squared Error (RMSE) Divided by 197 
the Mean of the Actual Values 

Mean Absolute Percentage Error (MAPE) Bie 


(in percent) 


Actual Values 56128 59995 59916 59795 
Predicted Values 49737 49396 45437 46556 
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the OLS model indicated the presence of only AR(1l) in tne 
residuals. The data was transformed and the model 
reestimated, resulting in a model with no indication of 
autocorrelation nor any violation of the assumptions. The 
results of the final model were a large improvement over the 
initial one but the predictions were once again poor. The 
predictions were actually better than expected because the 
cost center has at the most only one identifiable variable 
indirect position out of approximately 120 total indirect 
personnel in that department. The majority of the indirect 
labor would by definition not be expected to vary with the 
changes in the direct labor. 
4. Cost Center 300 

Table 8 shows the results for the NAVAIR Engineering 
Support Office (300). Very poor results were obtained for 
the initial model. The data was transformed to eliminate 
the presence of the AR(1) process and the model was 
reestimated. It was necesSary to transform the data again 
to eliminate the presence of the AR(4) process from the 
residuals. Analysis showed the residuals finally appeared 
to be free of any autoregressive process and none of the 
assumptions seemed to be violated. The results of the final 
model are quite good and should provide reliable 


predictions, 
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TABLE 8 
RESULTS FOR COST CENTER 300 


Model: ID300 = a + b DIRECT 


Untransformed Data 


standard Error of the Regression: 397 
Adjusted R-Squared: iO 
F-statistic (degrees of freedom): foi G30) 
Estimate of a: 26415 
Seanaard Error: 9615. 
T-statistic Ke ol 
Estimate of b: 02 
Standard Error: rg 0 
T-statistic PARAS 
Durbin-Watson Test Statistic: 63 
Niece eStmotatas tic: e oo 
Estimator for First Transformation (A ) 6870 
Estimator for Second Transformation (A, ) ~6596 


Transformed Data 


Standard Error of the Regression: 4071. 
Adjusted R-squared: SO 
B=Statastic: 63.6 
Estimate of as: 564. 
Standard Error: ge2 
T-statistic o4 
Estimate of b: Se 
Standard Error: sUl0'S 
T-statistic i296 


Prediction Results 


Correlation Coefficient (CC) Between -.626 
Actual and Predicted Values 

Root Mean Squared Error (RNSE) Divided by 307 
the Mean of the Actual Values 

Mean Absolute Percentage Error (NAPE) 4.7 


(in percent) 


ACeviaw Values S60 0mmeo col Gee 3/357 33103 
Predicted Values Biba 2 eo 956035731395 36.606 
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5. Cost Center 300 
The sequence followed for the model of the 
Production Planning and Control Department (500) was the 
sane as for the 300 department. The data was adjusted for 
the presence of AR(1) First and then for Ak(4) resulting in 
residuals that appeared to be free from the presence of 
autoregression nor violate any of the necessary 
assumptions. The results, shown in Table 9, are very good 
and should also provide good predictions for the uSer. 
6. Cost Center 609 
The results for the Production Engineering 
Department (600) are presented in Table 10. The initial 
regression produced very poor results. The data was 
transformed for the AR(1l) and then the AR(4) process and the 
model reestimated to give a fairly good resulting model that 
did not display any autoregression present in the 
residuals. The final model met all the required 
assumptions. The predictive ability displayed by the final 
model is actually very good as over ninety percent of the 
600 cost center's indirect is fixed indirect which does 
would not vary with changes in the direct workload. 
7. Cost_Center 620 
The results for the Plant Services Division (650) are 
shown in Table 1l. The results from the OLS were fairly 
good and these results could have been misleading due to the 


presence of autocorrelation. The data was transformed for 
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TABEE 9 
RESULTS 


MoGgel:s ID500 = 


Standard Error of the Regression: 
Adjusted R-squared: 
F-statistic (degrees of freedom): 
Estimate o€ a: 

Standard Error: 

Te-statistic 
Estimate Of Db: 

Standard Error: 

T-statistic 
Durbin-“atson Test Statistic: 
walters Test Statistic: 
Estimator for First Transformation 
Estimator for Second Transformation 


Standard Error of the Regression; 
Adjusted R-squared: 
F-statistic: 
Estimate of a: 
Stan@aca Error: 
T-statistic 
Estimate of b: 
Standard Error: 
T-statistic 


Correlation Coefficient (CC) Between 
Actual and Predicted Values 

Root Mean Squared Error 
tne Mean of the Actual Values 


tlean Absolute Percentage Error (NAPE) 
(in percent) 

Actual Values 210866 193054 

Predicted Values 2760 1 oeezl S07 


69 


(RMSE) Divided by 


PORMCOUST Centr enes00 


a + b DIRECT 


Untransformed Data 
40076. 
soe 
iG 6 
=o ao. 
Giooce 
= 3 
.5 
-08 
ao 
-68 
86 
(A ) .662 
(3) aa 


(1,32) 


Transformed Data 
24844. 
~48 
6) rae 
18043. 
72345 
2.49 
197 
03> 
SOS 


Prediction Results 


2930 
07 


6.6 


182240 
195081 


193472 
mO2295 


TABLE 10 


RESULTS FOR COST CENTER 600 


(lodel: ID600 = 


get D DIRECE 


Untransformed Data 


Standard Error of tne Regression: ed Clie, 
Adjusted R-sSquared: ~09 
F-statistic (degrees of freedom): 2g (1,28) 
Estimate of a: 67378 
Standard Error: 30829. 
T-statistic 2.0 
Estimate of b: ~064 
Standard Error: o03)7 
T-statistic liek 
DurdDin-Watson Test Statistic: 453 
Wallis Test Statistic: . a a) 
Estimator for First Transformation (A, ) TT 
Estimator for Second Transformation (A, ) 23770 
Transformed Data 
Standard Error of the Regression: BOs 21 
Adjusted R-squared: - 36 
F-statistic: 15.6 
Estimate of a: 9161. 
Standard Error: 3043. 
T-statistic $201 
Estimate of b: ~070 
Standard Error: -018 
T-statistic 3.95 
Prediction Results 
Correlation Coefficient (CC) Between ~.436 
Actual and Predicted Values 
Root Mean Squared Error (RMSE) Divided by a 
the Mean of the Actual Values 
Mean Absolute Percentage Error (‘IAPE) Ba 
(in percent) 
Actual Values 123.800 @L16589 123) 500 Zo 555 
Predicted Values 122480 117170 103918 105348 
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TABLE ll 
RESUBTSTSOR PeOST CENTER” 050 


Model: ID650 = a + b DIRECT 


Untransformed Data 


Standard Error of the Regression: 15644. 
Adjusted R-squared: 0 
F-statistic (degrees of freedom): 20.1 Cr, 2) 
Estimate of a: ~3696. 
Standard Error: 26644. 
T-statistic -.14 
Estimate of b: 146 
Standard Error: 2032 
T-statistic ao 7 
Durbin-Watson Test Statistic: ~41 
Wallis Test Statistic: 259 
Estimator for the Transformation (2) 793k 


Transformed Data 


Standard Error of the Regression: BES) 
Adjusted R-squared: 2oo 
F-statistic: 7 ave 
Estimate of a: Do fee 
Standard Error: 227205 
T-statistic 2.69 
Estimate of b: ~099 
Standard Error: ZOOL? 
T-statistic Bo > 


Prediction Results 


Correlation Coefficient (CC) Between 20 
Actual and Predicted Values 

Root MNean Squared Error (RNSE) Divided by ~04 
the Mean of the Actual Values 

Mean Absolute Percentage Error (iiAPE) a 


(in percent) 


Actual Values 104160 103673 98220 93905 
Predicted Values HealsZ2—) 103901] 97176 102108 
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the AR(1l) process and the model reestimated. The final 
model's residuals showed no presence of autocorrelation and 
the necessary asSumptions appeared to hold. The prediction 
results from the final regression model were excellent. The 
fi~odel should provide very good predictions into the future. 
8. Cost Center 100 

Table 12 presents the results for the Material 
Management Department (700). Analysis of the residuals of 
the initial regression of the untransformed data indicated 
that the presence of AR(1). The data was transformed and 
the model reestimated. The resulting GLS solution yielded 
residuals that showed no presence of any autocorrelation and 
no violation of any of the required assumptions. Although 
no predictive analysis was performed, all the statistics 
indicate that the final model would produce very good 
predictions. 

9. Cost Center 330 

The results for the Metal and Process Division (930) 
of the production department are provided in Table 13. 
Fairly good results were obtained for the initial regression 
which could have been misleading if the presence of 
autocorrelation had not been looked for. The residuals 
showed that they contained AR(1) so the data was 
transformed. The GLS regression was the only reestimation 
required to obtain residuals that did not indicate the 


presence of autocorrelation and none of the required 
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RESULTS FOR COST CENTER 700 


ttodel: ID700 = a + b DIRECT 


Standard Error of the Regression: 


Adjusted R-squared: 


F-statistic (degrees of freedom): 


Estimate of a: 
Grandarcd Error: 
T-statistic 
Estimate o€ b: 
Steamagard Error: 
T-statistic 
Durbin-Watson Test Statistic: 
Wallis Test Statistic: 
Estimator for the Transformation 


Standard Error of the Regression: 


Adjusted R-squared: 
F-statistic: 
Estimate o€ as: 
Standara Error: 
T-statistic 
Estimate of b: 
Standard Error: 
T-statistic 
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Gntransformed Data 
3889. 
23 
4.65 
30365. 
20617. 
i, 47 
056 
O26 
Re: 
8758 
oe 
(P, ) oad 


ips.) 


Transformed Data 
24235 
~ 89 
96.1 
10691. 
25104. 
4.14 
066 
0067 
9.8 


TABLE 13 
RESUETS FOR COS PeeENTER@ oe 


Model: ID9307= a + beni kecr 


Untransformed Data 


Standard Error of the Regression: 7766. 
Adjusted R-sSquared: soe 
F-statistic (degrees of freedom): 29 ieee (l255 
Estiinate of a: -24095. 
Standard Error: 16051. 
T-statistic “1.5 
EStiniaGe Ob sor ~104 
StanGdamdy Error: 019 
T-statistic 5. 4 
Durbin-‘iatson Test Statistic: ~ 84 
Wallis Test Statistic: 1.17 
Estimator for the Transformation (2 ) 25789 


Transformed Data 


Standard Error of the Regression: 230. 
Adjusted R-squared: 25 
F-statistic: SiO ke. 7 
Estimate of a: 4839. 
Standdead enor: 4106. 
T-statistic 1.18 
Estimate of b: O00. 
Standdtaseccor. moll 
T-statistic S50 


Prediction Results 


Correlation Coefficient (CC) Between ~ 15 

| Actual and Predicted Values 

Root tean Squared Error (RISE) Divided by ~066 
the Mean of the Actual Values 

Mean Absolute Percentage Error (lHAPE) ie 


(in percent) 


Actual Values 54044 47995 51345 48517 
Predicted Values 53154 53261 47348 48373 
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assumptions seemed to be violated. The prediction 
capabilities of the final model are very good. 
10. Cost Center 940 

Table 14 shows the results for the Avionics 
Division (940) of the production department. The results of 
the OLS regression were poor. The R-squared statistic was 
very low and the standard errors of the regression 
coefficients were very high. The residuals were analyzed 
and showed the Durbin-ijatson statistic to be significant, 
indicating the presence of AR(l). The data were transformed 
and the model reestimated. Analysis of the residuals showed 
no presence of autocorrelation nor violaticn of the required 
asSumptions. The predictive analysis showed that the final 
model 1S not as good as most of the previous ones but the 
predictions were fairly good. The model should produce 
reliaole predictions. 

ll. Cost Center. 920 

The results for the production department's 
Airframes Division (950) are presented in Table 15. The 
sequence and explanation of results are the same as for cost 
center 940. The residuals showed no presence of 
autocorrelation and none of the required assumptions 
appeared to be violated. This model can ve expected to 


yield fairly reliable prediction results. 
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TABLE 14 
RESULTS FOR COST MGR TER a2 0 


fodel: ID940 = a + b DIRECT 


Untransformed Data 


Standard Error of the Regression: locus. 
Adjusted R~Scquared: 2D 
F-statistic (degrees of freedom): LOS (azo 
Estimate o€ a: 7097. 
Standard Error: 29618. 
T-statistic ~24 
Estimate of b: oleae) 
Standard Error: ~036 
T-statistic B25 
Durbin-Watson Test Statistic: ~ 49 
Wallis Test Statistic: nO 
Estimator for the Transforination (2 ) ~Po03 


Transformed Data 


Standard Error of the Regression: 10974, 
Adjusted R-Squared: 044 
F-statistic: 23.6 
Estimate of a: 6945. 
Standard Error: 4500. 
T~statistic deg. 
Estiimate of b: ~090 
Standard @Eenor ;: ~019 
T-statistic 4.86 


Prediction Results 


Correlation Coefficient (CC) Between 64 
Actual and Predicted Values 

Root Mean Squared Error (RMSE) Divided by o13 
the Hean of the Actual Values 

Mean Absolute Percentage Error (MAPE) 12h 


(in percent) 


Actual Values 95903 75818 69715 72570 
Predicted Values 86819 87959 79602 81190 
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TABGE 15 
REGUGTo Por COSl cbiaer J950 


N@Gel. 3bos0 = a7 bp DIRECT 


GUntransf£orined Data 


Standard Error of the Regression: oO 
Adjusted R-Squared: eo 
F-statistic (degrees of freedom): a5 (is 2.) 
Estimate of a: =——94 057, 
Standard Error: 30445. 
T-statistic as Lr 
Estimate of b: see 
Standard Error: ~036 
T-statistic Beek 
Durbin-Watson Test Statistic: 045 
Weamelis Test Statistic: Loe 
Estimator for the Transformation ee 77 SY 


Transformed Data 


Standard Error of the Regression: EOS. 
Adjusted R-sSquared: 2OO 
P=Statistic: L726 
Estimate of a: 745. 
Standard Error; SOs 
T-statistic ate 
Estimate o£ b: 068 
Standara Error: OMG 
T= Statistic 4.2 


Prediction Results 


Correlation Coefficient (CC) Between oS 
Actual and Predicted Values 

Root Mean Squared Error (RMSE) Divided by ~14 
the Mean of the Actual Values 

flean Absolute Percentage Error (liAPE) e2ee3 


(in percent) 


Pecudal Values 63186 72649 60426 48589 
Predicted Values Doge —o9goGa. of583 54756 
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12. Cost Center 260 

Table 16 contains the results of the regressions 
for production's Power Plant Division (960). The initial 
regression yielded fairly good results which, again, could 
have been misleading if the residuals had not been examined 
for the presence of autocorrelation. ReSiaual analysis 
indicated the presence of both AR(1l)and AR(4). The data was 
first transforined to eliminate the presence of the AR(1l), 
and then the model was reestimated. Analysis of the 
residuals From the GLS regression showed ne indication of 
the presence of autocorrelation nor any viclation of the 
required assumptions. The model can be expected to give 


excellent predictions into the future. 


Di oO Ae 

Regression models were obtained for each of thirteen 
NARF cost centers to help management predict their required 
indirect for four quarters into the future conditioned on 
the direct workload for the respective quarters. The 
analysis resulted in obtaining reliable, useful inodels for 
all but the 00/01 cost center which does not have any 
variable indirect in it, only fixed indirect. Table 17 
contains a summary of all twelve of the models. 

The first column of the table presents the percent of 
total NARF indirect cost each cost center modeled contains. 


The twelve models presented provide reliable predictions for 
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ZABLE 16 


RES UBEomeORsGOsat GCEYNEER 960 


[lodel: ID960 = a + b DIRECT 


Standard Error of the Regression: 
Adjusted R-squared: 
P-statistic (degrees of freedom): 
EStimaiwe of a: 

Standame Error: 

T= Giese 1S take 
Eotummacte Of b: 

Standard Error: 

T-statistic 
Perrolifewatson Test Statistic: 
ewer s Test Statistic: 
Eocimaeor £O0r the Transformation oom) 


Stendard Error of the Regression: 
Adjusted R-Ssquared:; 
P=staitiscitkac: 
Estimate o€ a: 
StanG@ard Error: 
M~statistic 
EStimate of b: 
Seandard Error: 
T=staeiseac 


Correlation Coefficient (CC) Between 
Actual and Predicted Values 


POoe ean Squared Error (RHSE) Divided py 


the lIean of the Actual Values 
Mean Absolute Percentage Error (NAPE) 
(in percent) 


Boo 
33.223 


Bee 
34548 


BRetual Values 
Predicted Values 


19 


30062 
oO 


Untrcansformed Data 
Gio 7 4a 
~ 46 
Dee) 
=20229. 
ih2 265. 
= les. 
Or 
~015 
Saw 
On 
IG 
eo s4 


(1,28) 


Transforined Data 
4807. 
043 
252 
2502. 
3419. 
eli 
.0d3 
~009 
4.8 


Prediction Results 
~66 
ore 


Zo 


32611 
32520 
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SUMMARY CONTAINING FINAL WODELS 


All the final models are presented below. 


The variable 


DIRECT 1s the summation of production direct hours (930 + 
940 + 950 + 960) in a quarter. These models all use data 


that has been transformed as described herein, 


Final Percent Adjusted CE 
llodel Oe R-squared 

Indinect 

Cost,FY871 
1D012=12823+.028DIRECT ee ool = 016 
ID200=10495+.043DIRECT Grr ~64 =U 
ID300= ~O41DIRECT Seu On =e Gs 
ID400= ~O87DIRECT Ta ~ 82 024 
ID500=18043+.197DIRECT Tc ~48 a3 
ID600= 9161+.073DIRECT hoe 0 ~.43 
ID650= 5979+.099DIRECT TO 69 OU 
ID700=10691+.066DIRECT Tras 89 NA 
ID930= ~O61DIRECT 4.7 oe le 
ID940= ~O9O0DIRECT 6.6 ~44 064 
ID950= ~0638DIRECT 4.9 oo 73 
ID960= ~O43DIRECT 35 » 43 66 
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RMSE 
NEAN 


04 


MAPE 
(%) 


hee 
tN 0 0) & Boh 


—~e © @ @ e *® @ 
WO WWND YD UI WO) = ~) CWO 


over 92 percent of the total NARF indirect labor. In 
addition to the final models, the adjusted R-Squared, the 
correlation coefficient (CC) between the actual and 
predicted values, the root mean squared error (RMSE) divided 
by the mean of the actual values, and the mean absolute 
percentage error (NAPE) are also Summarized in the table. 
The intercept term in each equation represents the fixed 
indirect within the cost center being modeled. The absence 
of an intercept term implies that cost center has little 
fixed indirect in comparison to the amount of variable 
Indirect. Thus, the misSing intercept term in the 930, 940, 
950, and 960 production cost centers was not unusual. Cost 
centers 300 and 400 do not contain an intercept term in 
their models either. The data used in the analysis were 
adjusted over time to eliminate the effects of 
Organizational changes. These adjustments account for 300 
and 400 cost center's models not having an intercept term, 
as they were both significantly affected over time. Both 
the 300 and 400 models can be expected to provide excellent 
predictions as previously mentioned and presented in Table 


By) 
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Viz CONCLUSZONS 


The objective of this thesis was to collect and analyze 
data from the Naval Air Rework Facility (NARF) in Alameda to 
develop aS many as seventeen forecasting models. Each model 
would pertain to a specific cost center within the NARF and 
be used by management in their deciSion making process to 
forecast future indirect labor requirements. The models had 
to be understandable, easy to use by personnel from varied 
backgrounds, and reliable. These objectives required the 
study of the NARF organization, their budget planning 
procedures and constraints, their accounting methods, and 
their statistical recording procedures. The objectives also 
required a detailed study of autocorrelation; its effects 
on regression residuals, OLS regression results, and the 
elimination of the autocorrelation process that is present 
in the residuals. The final objective was to conduct 
predictive analysis on the final models that were presented 
in order to evaluate their forecasting reliability. 

It was determined that four cost centers did not need a 
forecasting model (02/800, 900, 903, 050) due to the small 
number of indirect personnel in each of those cost centers, 
and the small fluctuation of indirect worked in each 


quarter. It was not possible to obtain a reliable model for 
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Eaemu O07 0) Cost center. Pht Beet 1S ne tedisturbing 
because the cost center has only fixed indirect personnel 
whieh by déetinition do net wanyewith the direct workload. 
The summation of production's direct labor hours in each 
quarter was used as the independent variable for all models. 
Of the remaining major cost centers it was determined 
that the final models for 012, 200, 940, and 950 should 
yield fairly good forecasts into the future. All the 
remaining cost centers' models should perform very well in 
assisting management in their decisions for indirect 
personnel requirements. It is imperative that the user 
understand that the models presented use data that has been 
transformed to eliminate the effects of autocorrelation in 
the residuals of the regression. The values of the 
estimators of the correlation coefficients used for the 
transformations were provided with the presentation of the 
results for each model. The value that is forecast is 
indirect hours that are worked. Expected leave and time 
allowed must be added to the predictions to obtain a figure 
which equates to total personnel required in that quarter. 
The analysis of the predictions for leave and time 
allowed is an excellent area for thesis research. A study 
of the reliability of the CWPABS program is also warranted. 
There are currently large variances between the predictions 
the model is making for direct workloads and the work that 


is actually being performed. It was noted that the cost 
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center statement (7310-68) is no longer provided to the NARF 
in microfiche form. The statement provides valuable 
statistical data that should be available to the 500 
department. The policy of retaining that data on microfiche 


should be reinstated. 
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APPENDIS 8B 


COMPUTER PROGRAMS 


Diora pencixeCOncalns a wlrsting Of the programs that 
were used in the analyses performed. All of the functions 
are written in APL and contain documentation. 

Tne initial regressions were performed using the OLS 
function. Next the GLS fEFunction was used to perform a 
regression on data tnat had been transformed for either the 
AR(1) or AR(4) process uSing the TRANS Function. Both the 
OLS and the GLS functions are based on a function written by 
Husygrave ang Ramsey (1981, ppy.254-258). The predictive 
analysis was accomplished by uSing the PREANAL cEunction 
unless the Final model included the intercept term, in which 
case the PREANALINT function was used. The final 
computations and output format for the predictive analysis 


were performed with the PRED function. 
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APPEUDIZ 


DATA 


This appendix contains the data for the major NARF cost 
centers. An entry of -99 indicateS there was no obServation 
listed on the microfiche for that cost center during that 
quarter. Data are listed by direct hours worked for all 
cost centers, followed by adjusted indirect hours worked, 


and leave hours charged. 
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157457 
148371 
134997 
1600 
170745 
171481 
Po e2 
149315 
155389 
Oo Zoro 
13142109 
137464 


940 


221062 
247472 
25107 
227692 
200430 
249980 
239284 
218471 


200467. 


224529 
219038 
232576 
198632 
235500 
226580 
224866 
200666 
237 Jou 
228307 
231303 
192323 
244240 
232044 
222084 
196808 
220162 
227247 
218222 
194662 
204946 
201232 
218460 
192577 
197376 


950 


206739 
243019 
241237 
211656 
179589 
242255 

70088 
230257 
Ze 317 
263848 
248494 
293900 
243190 
ga 1179 
277563 
257268 
212831 
266391 
255524 
228707 
194788 
218137 
217361 
238519 
225249 
250543 
243238 
249571 
213473 
221059 
220081 
219219 
213232 
226366 


960 


dA 7 
216554 
2aS2 2 
185388 
143622 
220770 
198289 
189027 
177696 
199871 
195406 
232008 
193728 
Zo 413 
232897 
233397 
180762 
193308 
185539 
192075 
160100 
196013 
192269 
187372 
179847 
207162 
220956 
229600 
190871 
198475 
203642 
229012 
Z016)/ D9 
214005 


ANY 


050 


148 


QDOODOODOONTDWAOAOCMOCIOAOIGAOIGOINIAIDAIDAIAANAANAAOAAOOOO FOOWOAOO 


02/800 


1499 
2395 
3042 
2969 
1192 
sii3 
3852 
3145 
3164 
3340 
doo 
3306 
287-5 
2818 
2078 
2854 
2623 
2589 
3319 
3528 
9612 
6407 
3937 
9325 
4895 
5322 
UQUIL 
9895 
5606 
6255 
6006 
6305 
9772 
6619 


154 


Darect. Hours: 


Bk 


791 
792 
733 
794 
801 
802 
803 
804 
ead 
612 
813 
814 
8121 
S22 
823 
824 
831 
832 
833 
834 
Sui 
842 
843 
844 
Se 1 


Siac 
854 
861 
862 
863 
864 
871 
872 


012 


00/01 


~ 99 
~ 99 
~ 99g 
66 


200 


Sai 
5999 
6257 
6675 
6203 
22-3 
S046 
2746 
2400 
2893 
2628 
4223 
4633 
3536 
Joero 
3875 
4160 
322 
9789 
5486 
6714 
6480 
6729 
oe 
6279 
8476 
8020 
9415 
7490 
9988 
11133 
10209 
9086 
6458 


300 


37.045 
41521 
45370 
U564U8 
e4025 
65659 
65469 
69824 
55944 
70805 
74286 
77932 
60993 
LOW 7 
69799 
72149 
a2 2a 
66983 
67050 
66296 
36662 
73220 
T2028 
70280 
62747 
69722 
Viele 
69810 
39622 
65493 
67241 
75699 
68250 
73554 


400 


A 2Gi2 
1515 
1576 
2200 
1149 
1421 
1088 
705 
USYy 
372 
326 
“O06 
ae 
531 
566 
356 
140 
u80 
385 
326 
152 
343 
u32 
138 

56 
5381 
UUT7 
546 
363 
22u 
aon 
769 
605 
338 


3 


500 


41872 
50603 
47563 
37208 
33880 
46280 
49297 
92280 
U74U77 
62611 
58402 
99297 
4U4U819 
2253 
394632 
93817 
46289 
a0 
94376 
31905 
48670 
62383 
68777 
PD aa 
91919 
99676 
a7 7 37 
60897 
49708 
97914 
69687 
77958 
69566 
736 


600 


91235 
9847 
8128 
12033 
7789 
4875 
1914 
6540 
U5S4Q 
7889 
169 37 
Li 7 8 
9286 
ato 12 
16788 
11724 
10382 
14991 
aero, 
19644 
13432 
13450 
13322 
13989 
1253 
13810 
16786 
=) Zon 
16949 
17038 
19142 
19879 
Pores 
19865 


650 


14736 
13791 
15257 
12375 
13406 
16387 
14699 
11910 
12734 
6176 
4078 
2740 
2188 
U320 
1619 
1681 
760 
1979 
1246 
662 
a2 SS 
99u 
2144 
4351 
1623 
1886 
3072 
u53 
1348 
2250 
2227 
1929 
2642 
3739 


700 


Ga 


NO 
NO 
OOwWoNnNWOWDO00RD00 0 


G) 


Gd) 
i= 
~J 
© 


6197 
248 
90 


UGU 


608 
Oe 
1932 
1470 
1380 
1556 
370 
1481 

999 

746 
1389 


Indirect Hours: 


YR 


791 
792 
793 
794 
801 
802 
803 
804 
811 
Sk 
813 
814 
821 
822 
823 
824 
Sst 
832 
833 
834 
841 
842 
843 
BUY 
851 
852 
853 
854 
861 
862 
863 
864 
S77 
872 


930 


68360 
62601 
91962 
80025 
69063 
88075 
69528 
54025 
99931 
65253 
66763 
69595 
68342 
84757 
87 350 
83073 
67133 
7 te 6 
72444 
70651 
58477 
65193 
94706 
93894 
49736 
92002 
96027 
60420 
UQU82 
395906 
SHOWY 
47995 
51345 
48517 


940 


111437 
140020 
138900 
115395 
107892 
126066 
116919 
de 1O7 
102865 
125750 
109673 
114523 
An 2 Sie 
P20 
io 129 
128962 
114704 
125644 
117580 
106253 
105354 
116580 
81370 
86992 
82266 
92259 
89089 
86458 
78594 
88875 
95903 
75818 
69715 
72370 


950 


36510 
47821 
U7964 
35929 
34885 
49122 
8495 
48792 
49508 
Joo 16 
58437 
73071 
78335 
90761 


109288 


97201 
79236 
80292 
73765 
694UUS5 
B22 7 
79050 
60400 
40677 
36208 
47479 
98154 
54766 
52475 
99992 
63186 
72649 
60426 
48589 


960 


42930 
54705 


54635 
41400 
33 7G 
48933 
4uO890 
38907 
39741 
UU733 
46390 
55595 
53326 
55718 
57431 
58266 
47163 
54399 
49154 
43670 
41014 
45508 
JO UGS 
28466 
32877 
37932 
41598 
32200 
30561 
33206 
83267 
32278 
30062 
32604 


114 


050 02/800 
422 ~ 9g 
U14 ~ QQ 
UU 1 ~ QQ 
U36 ~9Q 
402 64 
260 oz 
0 a2 

~ 9g Ug 
“gg 0 
“9g uO 
~ QQ uO 
“gg 78 
“gg 7ga 
~ QQ 1065 
~Qg 1667 
“gg 3465 
~ QQ 37a 
0 3926 
2645 933 
3599 501 
61 U24U 
S23 398 
5uu7 594 
31 25s 
3406 4u58 
2213 U55 
1672 1098 
2153 1185 
1580 735 
2061 720 
1799 964 
2214 850 
1904 760 
3399 872 


Indirect 


1g hs 


791 
792 
793 
794 
801 
802 
803 
804 
i 
812 
813 
814 
821 
822 
823 
824 
831 
832 
833 
834 
841 
842 
843 
B44 
Si 1 
Si 2 
S53 
854 
861 
862 
863 
864 
a7 i 
872 


HOUrS: 

O12 00/01 
U7989 6263 
49210 8743 
49527 9391 
49266 8840 
49616 9715 
49926 10589 
50693 12376 
Shoo Ss 33 oS 
50u51 be2 1 Ay 
51162 £3713 
51943 i 5) aro. 
§2415 16374 
51405 14209 
52396 16178 
Sa6o0 4UU236 
65048 47832 
61284 39060 
62809 41938 
50776 23855 
55026 222617, 
Doo 19006 
54354 20632 
55032 22666 
54278 20920 
54U8uU7 22629 
54891 19600 
55989 26600 
58791 26146 
438326 17302 
48759 23164 
53500 23433 
49303 25140 
48932 I) lair) (3500 
54Qu49Q 27186 


200 


a2 781 
6307 2 
99658 
aco cm 
oii 
62876 
98502 
58415 
90819 
a) 62 
99820 
S20 7-4 
594818 
67566 
49196 
40199 
36592 
4UUu857 
99379 
361516 
47935 
94560 
92204 
96811 
45559 
90993 
92366 
50286 
48075 
54798 
J62 26 
99995 
99916 
99795 


300 


49135 
he? 
98379 
T3979 
a7117 
41908 
41091 
41854 
38849 
41533 
42233 
43680 
4OO4S 
41688 
41085 
41796 
36898 
36806 
396005 
39001 
31640 
34621 
33759 
34225 
30714 
est) yal 
36158 
37547 
38451 
37589 
36633 
32638 
ago 07 
38103 


eS) 


400 


70471 
e729 
80574 
77108 
69584 
713182 
67986 
67 5 oA 
67485 
84226 
86587 
86740 
71474 
86485 
85593 
80346 
71839 
78819 
76354 
70834 
o8isi 
64526 
65231 
64920 
20206 
97174 
64262 
61066 
97888 
61969 
70078 
69305 
62138 
63621 


500 


264594 
305367 
325044 
268951 
278124 
290024 
292224 
286370 
265029 
293279 
314407 
305286 
289133 
329976 
342548 
351396 
298619 
299228 
291398 
266439 
223297 
256194 
228435 
211432 
188748 
198667 
234272 
216018 
204771 
223677 
210866 
193054 
182240 
193472 


600 


72268 
88877 
98975 
90468 
81198 
Sie os 4 
92727 
ga735 
76621 
81301 
96268 
a ays aio 
1.07 330 
126427 
131436 
Toes 3 
118549 
119230 
i279 oe 1 
1344usYy 
114443 
121708 
114620 
dpe 2 
101881 
115304 
31567 
29622 
a Ibo 7 
124794 
123800 
116589 
ZU 
126965 


Indirect Hours: 


WR 


791 
792 
793 
794 
801 
SoZ 
803 
804 
811 
812 
813 
814 
821 
822 
823 
824 
831 
832 
833 
834 
841 
842 
843 
844 
Sra 
852 
853 
854 
861 
862 
863 
864 
871 
872 


650 


ieee 8 
150340 
153945 
136609 
129690 
143992 
124852 
128618 
117405 
125948 
132331 
138989 
126885 
145834 
143376 
140403 
116148 
114443 
108096 
112612 
95650 
107722 
98355 
90264 
86834 
94801 
99431 
97024 
93891 
105742 
104160 
103673 
98220 
93905 


a 


700 
33873 
a207 

~ 8 


Leave Hours: 


mR 


aa 
792 
733 
794 
801 
802 
803 
804 
Sad 
SZ 
813 
814 
eZ 
822 
823 
824 
831 
832 
833 
834 
84ul1 
842 
843 
844 
S52 
852 
853 
854 
861 
862 
863 
864 
871 
872 


550 


99551 
4U4U980 
63954 
66603 
95639 
4u998u4 
43858 
504390 
63270 
39076 
40706 
a aio 
71061 
48180 
47499 
596780 
67017 
37244 
40593 
46227 
62496 
S702 1 
36647 
38544 
36/65 
gia 
33693 
42124 
94399 
3610) 02 
28713 
42304 
52998 
37154 


940 


99999 
97986 
63997 
75988 
26 te3 
58688 
56409 
PAG deez 
90265 
SoZ 07 
50824 
69523 
g0191 
52185 
939333 
67813 
88419 
“7481 
48790 
63072 
91969 
50417 
94108 
63990 
87678 
44378 
44702 
58275 
82598 
45575 
4Q514 
2022 
72923 
50228 


950 


78054 
ake ga (0 al 
591949 
59482 
72270 
47023 
Bea at 
98832 
84317 
50037 
46772 
60629 
92028 
97092 
71695 
60242 
88176 
90992 
47860 
96281 
77481 
43425 
45192 
92383 
73224 
42998 
42823 
I3 216 
71374 
49067 
39421 
53459 
70700 
58065 


960 


78106 
42626 
453952 
51086 
70042 
51094 
37990 
§1428 
66734 
41659 
34674 
48697 
69104 
38897 
38665 
JS.) 
64681 
31810 
35045 
44071 
63466 
37799 
37704 
46102 
66652 
339197 
36407 
43862 
60823 
38921 
28947 
43669 
98166 
47458 


IIE 


050 02/800 
146 565 
gg 527 
78 150 
83 Tek 
126 604 

52 UY 3 

0 340 
~gQ 311 
~ QQ S32 
~ QQ 229 
~gQ Aiea 
“QQ 330 
“QQ 529 
~Qg 226 
~ QQ 182 
~ 9g 242 
“99 1040 
as 2 263 
104 233 
145 528 
366 1005 
653 615 
338 1020 
4U55 3239 
830 1563 
164 1058 

88 1466 
215 ‘1435 
374 1167 
116 1811 
15.3 940 
408 1050 
594 1286 
353 1476 


og 


679 
916 


Leave Houcs: 


oR 


791 
792 
793 
794 
801 
802 
803 
804 
811 
812 
aia 
814 
821 
Sie 
823 
824 
831 
832 
833 
834 
841 
842 
843 
BUY 
851 
6 2 
Sis 
854 
861 
862 
863 
864 
871 
872 


O12) 007 oH 
1267 2327 
Zoo 1386 
G73 12a7 
540 997 
Lito 2258 
6cany 1257 
601 1180 
824 1668 
1292, 2976 
741 1540 
419 729 
LOSBy 252 
12239 2355 
7605 Waid 
1272 2896 
1584 3689 
3324 7756 
1397 3223 
827 1770 
1758 3899 
3149 6869 
1848 4054 
2329S Sit 
735 1408 
3333 7378 
1509 3286 
9888 1176 
7962 3699 
12620 3o23 
7iae8)) 1385s 
9350 2006 
7834 4746 
11407 5778 
9611 5643 


200 


11376 
9356 
a2 
9109 

10693 
2 Zee 
7100 
9121 
Joes 
4652 
3861 
6560 

12928 
8337 
6540 
7467 

10803 
4270 
6 S25 
7720 

15093 
7008 
8409 
9807 

13730 
7599 
6975 
9505 

15333 
9113 
8661 

13450 

18159 

13669 


300 


26871 
113394 
14683 
19088 
26077 
14995 
14904 
18146 
26658 
14786 
12184 
18146 
25487 
13365 
18398 
ti 
25421 
11933 
13002 
16184 
25314 
12163 
13109 
17840 
26634 
11345 
12418 
19587 
24953 
14579 
12420 
17082 
23936 
17949 


400 


22796 
11784 
T2203 
16039 
21658 
15697 
11996 
16249 
20338 
10986 
10009 
13621 
18452 
10679 

9657 
13381 
18287 

Sod 

9539 
12821 
17769 
10324 

8717 
12460 
17419 

9938 

7970 
10539 
15954 

9397 

S615 
Piso 
16225 
10036 


jee 


500 


99999 
96622 
J Zoe? 
67990 
97711 
67769 
98612 
71209 
98103 
977393 
47298 
64337 
84254 
53144 
92 5100 
68876 
91297 
44335 
43149 
94216 
79003 
4O9u1 
47838 
56605 
Pazas 
40906 
41059 
56100 
72418 
91082 
37912 
91200 
66510 
49419 


600 


27220 
ime7 3 
13499 
16467 
24780 
15369 
12979 
15763 
22764 
11781 
13031 
20124 
27677 
15 4173'S 
15368 
21404 
31931 
125123 
12501 
21725 
32161 
15899 
18065 
23132 
33423 
16479 
16078 
25767 
34456 
20553 
15088 
20976 
33523 
23140 


650 


48131 
24350 
30487 
35482 
48567 
28849 
24128 
30575 
36065 
22973 
20695 
24086 
31658 
21971 
19814 
26333 
32636 
19664 
19612 
ZileoiZ 2 
26480 
17267 
Liao 7 
21114 
26407 
15204 
13423 
18647 
24035 
17458 
14564 
2221.5 
24384 
17794 


700 


Oo000CaOaAOaO aAOoOaaoOeo0o0o0 °o 


1843 

3737 

9672 
Toso 
11073 
12809 
14940 
22161 
13901 
12030 
15737 
19978 
13023 
11391 
13058 
20809 
14408 
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